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Abstract 

The paper establishes tight lower bound for effective conductivity tensor K, of two-dimensional three- 
phase conducting anisotropic composites and defines optimal microstructures. It is assumed that three 
materials are mixed with fixed volume fractions and that the conductivity of one of the materials is infinite. 
The bound expands the Hashin-Shtrikman and Translation bounds to multiphase structures, it is derived 
using a combination of Translation method and additional inequalities on the fields in the materials; similar 
technique was used by Nesi (1995) and Cherkaev (2009) for isotropic multiphase composites. This paper 
expands the bounds to the anisotropic composites with effective conductivity tensor K, . The lower bound 
of conductivity (G-closure) is a piece- wise analytic function of eigenvalues of K m , that depends only on 
conductivities of components and their volume fractions. Also, we find optimal microstructures that realize 
the bounds, developing the technique suggested earlier by Albin et al. (2007) and Cherkaev (2009). The 
optimal microstructures are laminates of some rank for all regions. The found structures match the bounds 
in all but one region of parameters; we discuss the reason for the gap and numerically estimate it. 

Keywords: multimaterial composites, optimal microstructures, bounds for effective properties, structural 
optimization, multiscale 

1 Introduction 

The problem The paper investigates the lower bound for effective conductivity and optimal micro- 
geometries of three-material composites (plane problem). We assume that two mixing isotropic materials 
have finite conductivities ki and fe, (0 < fci < fe) and the third one is a superconductor k$ = oo, the volume 
fractions mi > 0, ni2 > and rriz = 1 — mi — vri2 > of the materials are fixed. The conductivity of a 
composite is characterized by an anisotropic effective conductivity tensor K* that depends on the properties 
of mixed materials and their volume fractions, as well as on microstructures. We describe the bounds 
of G-closure (Lurie and Cherkaev, 1981) - the set of all effective properties of composites with arbitrary 
microstructure. The G-closure boundary depends only on k\, ki, mi, m,2. Optimal microstructures realize 
the bound if their effective conductivity lies at the G-closure boundary. 

We find the bound solving a variational problem of minimization of K, with respect to microstructures 
(Section[5]). Namely, we apply two orthogonal external fields of different magnitudes to a periodic composite 
and minimize the sum J of the corresponding energies of the composite, varying the microstructure occupying 
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the periodicity cell fl. The computed value of J allows for computation of the add outer bound of G-closure, 
as discussed in Sections 12.21 and 12.31 The matching microstructures (minimizing sequences) are found by a 
different technique that was introduces in (Albin et al, 2007a, Cherkaev, 2009) and is described in Sections 
13.31 and [6] by assumption, optimal structures are laminates of some rank. The effective properties of the 
structures form the inner bound of G-closure. When the outer and inner bounds coincide, they are exact 
and the G-closure is determined. We show that our bounds are exact in all domains of parameters but one. 
In the last domain, we estimate the gap between the outer and inner bounds. 

Remark 1.1 The complementary upper bound can be established by a solution of a dual problem, in which 
conductivity ki are replaced by resistivity pi = 1/ki. In the considered problem, one of the component is a 
superconductor (k-j, = oo) which makes the dual bound trivial - the effective resistivity can be arbitrary large, 
or K^ 1 > 0. The obtained results allows for the upper bound determination for the G-closure of materials 
with conductivities ki — < &2 < &3 < oo. 

Bounds The problem of exact bounds has a long history. It started with the bounds by Voigt (1928) 
and Reuss (1929), called also Wiener bounds or the arithmetic and harmonic mean bounds. The bounds 
are valid for all microstructures and become in a sense exact for laminates: One of the eigenvalues of K* 
of a laminate is equal to the harmonic mean of the mixed materials' conductivities, and the other one - 
to the arithmetic mean of them. The pioneering paper by Hashin and Shtrikman (1963) found the bounds 
and the matching structures for optimal isotropic two-component composites, and suggested bounds for 
multicomponent ones. The exact bounds and optimal structures of anisotropic two-material composites were 
found in earlier papers (Lurie and Cherkaev, 1982, 1986, Kohn and Strang, 1983, 1986, Tartar, 1985) using a 
version of the translation method (see its description in books (Cherkaev, 2000, Allaire, 2001, Milton, 2002, 
Dacorogna, 2008)). The method is equivalent to building the polyconvex envelope of a multiwell Lagrangian, 
as it was shown by Kohn and Strang (1983, 1986); the wells are the energy of the materials plus their cost 
(here, "cost" is the dual variable to the volume fraction of material in the composite). The theory of bounds 
for the two-material composite is now well developed and applied to elastic, viscoelastic, and other linear 
materials, see for example the books (Lurie, 1993, Cherkaev, 2000, Allaire, 2001, Milton, 2002, Dacorogna, 
2008). 

Bounds for multicomponent composites turn out to be much more difficult. Milton (1981) showed that 
the Hashin-Shtrikman bound is not exact everywhere (it tends to an incorrect limit when mi — > 0), but 
is exact when mi is larger than a threshold, mi > g m . Milton and Kohn (1988) suggested an extension 
of the translation method to anisotropic multimaterial composites, computed the anisotropic bounds for 
multicomponent composites and the optimal structures. Nesi (1995) suggested a new tighter bound for 
isotropic multicomponent structures, and Cherkaev (2009) further improved it and found optimal structures. 
The method is based on the procedure suggested by Nesi (1995) that combines the translation method and 
additional inequality constraints (Alessandrini and Nesi, 2001). These two latest bounds coincide in the case 
k$ = oo that is considered here. 

This paper extends these bounds to anisotropic composites. As in the early paper by Kohn and Strang 
(1983), we investigate the case when one of the phase has infinite conductivity, which significantly simplifies 
the calculation. The method is based on constructing a lower bound for the composite energy accounting 
for Alessandrini and Nesi (2001) constraints. Because of the constraints, the translated energies-wells can 
become nonconvex but are still bounded from below, an improved bound corresponds to this case. The 
method is described in Section [3] the results are summarized in Section [4] The energy bound turns out to 
be a multi-faced piece-wise analytic function of the problem's parameters. Like the translation bound, it 
depends only on conductivities of the materials, their volume fractions, and the anisotropy of a homogeneous 
external loading. The energy bounds and related bounds for the G-closure are derived in Section [5] 

Optimal structures In the paper, we prove that multiscale laminates realize the G-closure bound. 
Similar structures - laminates of second rank - realize the G-closure bound for the two-material case (Lurie 
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and Cherkaev, 1982, 1986); three-material bound is achievable by more complex structures of the same kind. 
Optimal structures depend on the degree of anisotropy of the external loading. 

Remark 1.2 Generally, optimal structures are not necessary laminates: For example, Hashin and Shtrik- 
man (1963) first suggested "coated spheres" geometry, Milton (1981) and parallel coated spheres and later 
suggested (Milton, 2002) a method of transformation of optimal shapes, Lurie and Cherkaev (1986) sug- 
gested multilayer coated circles, Vigdergauz (1989), Grabovsky and Kohn (1995) and recently Lui (2008) 
suggested special convex oval-shaped inclusions, Gibiansky and Sigmund (2000) suggested "bulk blocks", Al- 
bin and Cherkaev (2006) proved the optimality of "haired spheres", and recent paper by Benveniste and 
Milton (2010) investigated "coated ellipsoids". All these structures admit separation of variables when ef- 
fective properties are computed. It is not clear yet if the laminate structure approximates any other optimal 
structure, see for example (Pedregal,1997, Briane and Nesi, 2004, Albin et al., 2007b). We show, however, 
that proper laminates are optimal for the considered problem. 

The topology of two-material optimal structures is simple and intuitively clear: for isotropic or moderately 
anisotropic loading, the less conducting material ki "wraps" the more conducting one k-2 (&2 > fci), so that 
&2 forms an nucleus and ki forms a core. If the anisotropy of the loading exceeds a threshold, the optimal 
structures degenerate into simple laminates. 

The multimaterial structures are more diverse and nonunique and require new ideas for constructing. 
Milton (1981), Lurie and Cherkaev (1985) and later Barbarosie (2001) described two types of isotropic 
structures that realize the multicomponent bound for sufficiently large volume fractions mi > g m of the 
weaker conductor ki < hi < . . . Later, Gibiansky and Sigmund (2000) expand the domain of applicability 
of Hashin-Shtrikman bounds to mi > g gs where the threshold g gs is smaller than the one previously known 
g g s < gm- They demonstrated new isotropic non-laminate microstructures (bulk structures) that realize 
this bound. Liu (2008) found another structures an optimal conductivity. Albin et al. (2007) extended the 
results of Gibiansky and Sigmund, (2000) finding anisotropic laminates that realize translation bounds for 
both isotropic and anisotropic structures in a range of parameters 

|fc*2 - k*l\ 

mi > g aC n, and - < g acn 

fc*l + fc,2 

where k,i and fc,2 are eigenvalues of K*. These inequalities restrict the range of volume fractions and degree 
of anisotropy of a composite that correspond to translation bounds. For isotropic composites (fc„i = fc„2), 
the range of applicability of the founded laminates coincides with the one of bulk structures founded by 
Gibiansky and Sigmund (2000). Structures that realize the isotropic bound for the whole range of volume 
fractions were found in (Cherkaev, 2009). 

In Section [6] we extent this result to anisotropic composites, finding new optimal structures that realize 
our new bounds. More exactly, we show that optimal laminates realize the bounds in all but one region. The 
topology of optimal structures depends on volume fractions of the mixing elements and loading anisotropy 
level. The structure adjusts itself to meet the sufficient optimality conditions that are found during derivation 
of the bounds. All the optimal microstructures are found by the same procedure suggested in (Albin et al, 
2007)) and based on (i) the energy bounds and sufficient optimality conditions for gradient fields inside each 
material, and (ii) the lamination technique that allows for satisfaction of these conditions. In all cases but 
one, the found laminate achieve the bounds, they are not unique. 

In the remaining case, the lower bound for G-closure is definitely not exact, hence the mentioned technique 
for building the structures is not applicable. In that case, we guess the best structures (that correspond to 
the upper bound of G-closure) basing on asymptotic behavior of optimal structures in neighboring regions 
and then numerically compute the gap between the structures and bound that is between the upper and 
lower bounds for G-closure). The gap is very small, see Section ^. 41 which shows that the suggested laminates 
(upper bound) and the lower bound accurately approximate G-closure. 
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2 The problem 



2.1 Equations and notations 

Consider a periodic composite formed by three materials. The materials ki occupy plane domains fii, i = 
1, 2, 3 C R.2 that form a unit periodicity cell fi 

(J fi< = fi, fi = {(11,12) : < si < 1, < x 2 < 1}. 

i=l,2,3 

The areas nij = ||fii|| of fii are fixed: mi +7712 +777,3 = 1, nii > 0, and there are no other constraints on fii. 

The conductivity of the composite is described by the system of equations that expresses the curlfrec 
nature of the field e, the divergencefree nature of the currents j, and the constitutive relation (Ohm's law) 
that joins the field and current 

V x e = 0, or e = — Vit, V • j = in fi, j = fae in fii (1) 

where u £ Hi (fi) is a scalar potential. Potential u and the normal current are continuous at the boundaries 
dik between fii and fi^, the conditions hold 

T-(ei-e k ) = 0, v-(ji-jk) = at d ik (2) 

where r and v are the tangent and normal vectors to dik- 

In order to determine effective properties of the composite, we subject it to two homogeneous orthogonal 
external loadings eo a and eob, 

lim e a (x) — e aa — ( „ ) and lim e b (x) = e ob = ( ° 
|z|->oo \0 J \x\^oo \r 

where < r < 1 and calculate the sum of energies. Thus, we deal with a couple of conductivity equations {TJ 
(denoted with subindexes a and 6) that differ only in boundary conditions. This couple can be conveniently 
viewed as a problem for a vector potential U = (u a ,Ub) T and 2x2 matrices E — (e a |ei>) and J = (j a \jb) in 
a cell fi: 

V • J = 0, E — -VU in fi, J = ki E in fii lim E{x) = E . (3) 

\x I — >oo 

E{x) is periodic in fi, J E{x) — Eq, Eq = ( J ^ (4) 

Energy In each material fcj, energy Wi — Wi a + Wo, of the coupled conductivity equation is defined as a 
sum of two energies Wi a and Wa, caused by the loadings eo a and eob, respectively, 

Wi = Wi a + W ib = ^hTr{EE T ) 

where Tr denotes the trace. The energy density of an anisotropic material with conductivity tensor K has 
a similar form 

W(K, E) = ~Tr(K E E T ) (5) 

Here, we assume that material 3 is a superconductor, £3 = 00, as in Kohn and Strang (1983). In fi3, field 
E is zero and the current is not defined. The energy of the superconductor is 
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The energy Wo(Eo) of the whole periodicity cell has the form 



1 2 
= ~Tr{K.E E%)= .^jgjE/ 



2 Easm(3),(4J JTJ in 

it defines the effective tensor A', and depends on parameters ki,mi,r and on subdivision fii 



Basis It is convenient to decompose the 2x2 matrix E using an orthonormal basis 

1 0\ 1/10 
a i = — f= I ^ , a 2 = 1 



2 V° 1/ \/2 V° - 1 

aa = 71 (-°i o) a4 = 7f(? J) (8) 



The decomposition takes the form 



where 



E = 7! Sl(jB)ai + 7f dl(£)a2 + 7! d2(i5)a3 + 7f S2(£)a4 (9) 



si(E) = ^=( £ ii + ^22), di{E) = ^=(Sn - £22) 

s 2 (E) = -^(E 12 -E 21 ), d 2 (E) = -L(E 12 + E 21 ). (10) 



The energy has a form 
Notice the representation for determinant 



W(k,E) = ^kTr{EE T ) = -*(«? + «2 + d? + <g). (11) 



det(E) = l(sl + s 2 2 -di-d 2 2 ). (12) 
In the next consideration, we assume that the average field is fixed: (si, s%, di, a^) € £ a u, where 



(s: /" si(x)dx = g i = ^"t-^ , / s 2 (x)dx = S 02 =0, 
I in v2 in 

f di(x)dx = D i = ^J^, [ d 2 (x)dx = D02 = oj (13) 
in v2 in J 



2.2 Bounds 



Energy Bounds We find a lower geometrically independent bound for Wo by arbitrarily varying sub- 
domains fii and but preserving their areas - fractions rrn of the materials in the composite. The bound 
depends on ratio r of the external loadings. Without lose of generality, we assume that the largest loading 
equals one, < r < 1. The case r = corresponds to only one loading, and r = 1 corresponds to two 
loadings of equal magnitude. The bound has the form, see 

B(E ki,mi)= inf W (Eo). (14) 

n^ : I n i =mi 

where Eo(r) is given by (|13|l . Below, we write B = B(r), assuming that other parameters are fixed. Notice 
that sum of the energies B stays constant if the order of external loadings is reversed, or labeling of the axes 
is reversed, or r changes its sign. 
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Algebraically, the bound B can be expressed through a tensor K* (Eo) of effective properties of the optimal 
composite with eigenvalues fc*i and k*2 and eigenvectors directed along OXi and OX2 axis, respectively. 



2B(E ,ki,mi,r) = k*i+k*2r (15) 



2.3 G-closure boundary 



The bound for the energy implies the bound for the effective properties tensor K t . To derive this bound, 
we use representation ([5]) 

2W comp (r) = (/e*i + k, 2 r 2 ) > 2B(r) (16) 

as a function of r. If the eigenvalues fc»i, fc,2 lie on the boundary of G-closure, the inequality for the bound 
becomes an equality, W CO m P (r) — B(r) = 0. 

Each value of r corresponds to an optimal pair (fc*i(r), fc*a(r)). To link these eigenvalues together, 
compute the enveloping curve 

[W comp (r) ~ B(r)] = 

of the r-dependent family of bounds and obtain the system for fc,i(r), k*2(r): 

, „ r dB , . 1 dB .„„, 

MO = 5-^, Mr) = j; F (17) 

Pairs fc*i(r), fc»2(r) form a parametric equation for the boundary of the G-closure set. If the parameter r 
can be explicitly excluded from system (|17[) . we obtain an explicit relation between r fc,i(r) and fc*2(r). 



3 Methods for bounds and optimal structures 
3.1 Formulation and procedure 

The technique for the bound derivation is described in (Nesi, 1995, ; Cherkaev, 2009). Here, we repeat the 
arguments of the last paper. As in the translation method (se, for example. Cherkaev, 2000, Milton, 2002) 
we construct lower bound using quasiaffmess of determinant (see Morrey, 1952 and Reshetnyak, 1967 ). 

/ det(VC/)dx = det(£ ) = 2r (18) 

JSi 

We rewrite {JJ adding and subtracting 2tdct_E = 2rt, where t is a real parameter, t £ R to the right-hand 
side terms 

Wo(E ) = inf [ (Tr(KEE T ) + tdet(E)) dx - 1 det(E ) (19) 

According to Translation method, we relax point-wise differential constraints on set of minimizers E = 
V(7, U £ Hi(n, R 2x2 ) + Eox in the integral in right-hand side of (|19[1 by replacing this set with a larger set 
£ of E: 

£ = {E: £eL 2 (fi,l 2x2 ), E &£ av ,} (20) 

We also use an extra constraint (Alessandrini and Nesi, 2001) det VU > a.e. in O, if det(£o) > 0, that 
further constraints the minimizers' set: E £ £+■ 

£+ = {E: Ee£, s\ + s\ - d\ - d\ > a.e. x £ Q.) , if r > 0. (21) 

The constraint is imposed as in (Nesi, 1995, Cherkaev, 2009). 
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Remark 3.1 The bound for isotropic composites obtained in (Cherkaev, 2009) uses a stronger inequality, 
which, however, coincides with $21]) for the considered here case hi, = oo. 

The relaxed problem for the lower bound does not contains differential constraints. It has the form 

Wo(E )>Y(r,t)-2rt (22) 

where 

Y(r,t) = inf / (k(sl + sl + dl+dl)-t(s\ + sl-dl-dl))dx (23) 

In the next section, we explicitly compute Y (r, i). Inequality (|22p gives a t dependent family of lower bounds 
on WX-Eo)- Let us call the best possible bound of this family B(r). We have 

W (r) > B(r), B(r) = max(F(r,t) - 2rt) (24) 

Using (|17p , we find the bounds for K,, or the bounds for G-closure. The constraint (|2ip could be slack or 
strict in different regions f2;, depending on the problem's parameters. If inequality (|21[) is slack everywhere, 
the procedure coincides with polyconvexification and gives the conventional translation bound. 



3.2 Calculation of the translated energy Y in the materials 

We split the integrals into a sum of integrals over f2; and minimize each term independently. This is possible 
because the differential constraint on the field is relaxed so that the boundary conditions between the fields 
in neighboring domains are neglected. Introduce the averages of the fields over Oi, 

Si, = — [ Sj (x)dx, Dij = — / d j (x)dx, j = 1,2. (25) 

and recall that the field in is zero: s 2 j(x) = d](x) = 0, Vi 6 SI3 => S 3J = D 3j = 0. Then (23J becomes 

Y{r,t)= min (V 1 + V 2 ), (26) 

O j j , iJl j to 

where 

S : {miS 1 !! + tt^&i = 5*01, m\Si2 + 1712822 = , 

miDu + m 2 D2i = D01, m 1 D 1 2 + m 2 D 2 2 = 0} (27) 

and 

Vi = inf / ((ki + t)s 2 + (hi - t)d 2 ) dx, i = 1,2, (28) 

s 2 >d 2 Jq. 

2 2,2 (J i2 1 i2 

s = S! + s 2 , d = di + d 2 
Here, Vi are functions of Sij , Dij , which we compute now 

Structure of minimizers Depending on the sign of ki—t, algebraic expressions for Vi and corresponding 
minimizers s(x),d(x) vary. We recognize three cases: 

(a) ki — t > is true. All terms in the integrand of Vi are convex. Applying Jensen's inequality to each 
of them, we have: 

V = mi(ki + t)(Sl + Sf 2 ) + rm(ki - t){D 2 a + D 2 2 ) (29) 
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The minimum is achieved when fields s, d are constant: 

sj(x) = , dj(x) = Dij a.e. x £ fij, j = 1,2 

(b) fci — t < is true. In this case, the d terms are concave, however the inequality holds (fcj — t)d 2 > 
(ki — t)s 2 . Replacement of the second term in the right-hand side in (|29p by (ki — t)s 2 decreases the integral 
value. We obtain: 

Vi = inf / 2k iS 2 dx = 2k i m i (S 2 1 + S 2 2 ) (30) 

where infimum becomes minimum again if minimizers s, d axe as: 

Sj = Sij, d 2 = s 2 a.e. x £ f2j or Sj = Sy, d 2 = S 2 X + S 2 2 = S 2 a.e. s £ fi 4 (31) 

The minimizers di(x) and d 2 (x) that correspond to the minimum assume the representations 

di(x) = Si cos(0(x)), d 2 (x) — Si sxn(9(x)), a.e. x £ f2j (32) 

where #(:r) is an arbitrary function. # does not affect the value of V (the bound (|30|) is independent of Di) 
but the averages Dn and Di2 depend on 9, they vary in a disk 

A 2 i + D 2 2 < (Si + S 2 2 ) . (33) 

All points in the disk correspond to some 8(x), the inequality becomes equality when 9 is constant a.e. 
x G Qi . 

Notice that 

det e = s 2 + s 2 — (d\ + d 2 ) = 0, a.e. x € fi, (34) 

In this case, constraint Q2ip play a crucial role because of the non-convexity of term (ki —t)d 2 . 

(c) If ki — t — is true, second term disappears, Vi becomes: 

Vi = inf / 2hs 2 = 2m i h(S 2 1 + S 2 2 ) (35) 
Jn 

which is similar to the previous case. However the minimum is achieved independently of d-fields that can 
vary in Qi keeping the relations d 2 < s 2 a.e.x £ fii, and s-components are constant, Sj = Sy, a.e.x £ fli. 

Remark 3.2 Accounting for inequality constraint, Vi can be written as 

Vi = $r i= ( inf s 2 id 2 g£ / n . ((h + t)s 2 + (h - t)d 2 ) dx if s 2 > d 2 
' 1 \ +oo otherwise 

The integrand of Vi(s,d) is a nonconvex saddle-type function of s and d. Due to the constraint, it is 
nonnegative and grows quadratically. One can check that expression for Vi in fS0\) is the proper part of the 
convex envelope of the integrand in Vi. The envelope is supported by the minimizers f31p . 

Summary The bound is derived in following steps. 

1. Arbitrary fixing translation parameter t £ [0, oo), conductivities ki and k 2 , and average field compo- 
nents Sij and Dij, find the lower bound (|28[) of translated energy Vi(t,k\,k 2 , Sij, Dij) in subregions 
Qi (see above). 

2. Fixing t and loading Eo(r), minimize the energy of the cell with respect to Sij, Dij and find optimal 
Sij (t, r, k\, k 2 , mi, m 2 ) and Dij(t, r, k\,k 2 ,rn\, m 2 ), and the bound Y(t, r, ki, k 2 ,m\,m 2 ) for translated 
energy. 

3. Find optimal t and energy bound B(r,ki, k 2 ,mi,m 2 ), using (|22ll . 

4. Express energy bound through effective tensor K, and find equation (|17|) for the G-closure boundary 
G(K,,ki, k 2 ,mi,m 2 ) = 0. 
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3.3 Optimal fields and optimal laminates 

Method for finding matching structures The exactness of obtained bounds is checked by finding 
a matching laminate micro-geometries that realize them; this is done regularly as in (Albin et al, 2007a, 
Cherkaev, 2009). In this section, we explain the method and optimal three- material structures. The method 
assumes that obtained bounds are exact so that the range of minimizers-fields inside each material is known. 
The exactness of the bound is demonstrated by constructing a sequence of multi-rank laminates, such that 
the field inside each phase strongly tend to fields in optimal ranges; consequently, the energy converges to 
the bounds, see for example (Conti et al,, 2003, Albin et al., 2007b) In the derivation of the bound, the 
ranges of optimal fields are not restricted by compatibility conditions ([2}, but are determined solemnly from 
the optimality of a relaxed problem (|22[) . The optimal laminate structure is a rank-one path between the 
optimal fields computed from sufficient conditions. If such path exists, the bound is realizable, and therefore 
optimal. 

A laminate with tangent r can connect two fields e\ and e 2 if they are rank-one compatible or det(ei— e 2 ) = 
0. In particular, laminates with tangent r = (1, 0) can connect two fields 

e 1 =f a n 1 ?) ande 2 =f a n 2 ^eigenvector 



bij \ b 2 

if a,2 — a\. The average field in the laminate is 



c6i + (l-c)6 : 



where c £ [0, 1] is the volume fraction of the first material. Here, we consider only orthogonal laminates, 
with r codirected with one of the eigenvector of e. The eigenvectors of laminates are the same. In the 
eigenvalues plane, each field is represented by a point vector (a, b). The compatible fields e\ = (ai,6i) and 
e 2 = (12,62) have a common first coordinate a\ — a 2 , and the set of laminates with all volume fractions 
is represented by an vertical interval ei2 = [ai,c&i + (1 — 0)62], c G [0, 1]. Similarly, the layers with the 
tangent (0, 1) connect fields with fields e\ = (ai, b) and e 2 = (a 2 , b), and the set of laminates with all volume 
fractions is represented by a horizontal interval [cai + (1 — c)a,2, b], c £ [0, 1]. 



Example 1. Second rank laminate (T-structure) The second-rank laminate can connect three 
fields 

ei = (ai,6i), e 2 = (a 2 ,b 2 ), e 3 = (a 3 ,&3), if ai = 02, 63 €[61,62] 
First, we connect fields ei and e 2 by laminates as in (|37[1 and choose c = Co so that 

Co 61 + (1 - c )b 2 = b 3 , 

then laminate ei2, (|37|) becomes compatible with ez and can be joined by an orthogonal laminate with 
tangent = (0, 1). We call the resulting structure 1/(12, 3) pointing out that fields ei and e 2 are joined by 
laminates, then the structure is rotated 90° (the rotation is shown by comma) and laminated with ez. This 
rank-one connection of ez with the laminate from ei and e 2 is possible if the relative volume fractions in 
that laminate are specially fixed. The volume fraction of the third material can be arbitrary. The average 
field ei2,3 in the structure is simply a sum ei2,3 = m\e\ + m 2 e 2 + mzez- 

Remark 3.3 Strictly speaking, the second-rank laminate is a sequence of structures, the corresponding ener- 
gies converge to the calculated values when the ratio of laminates scales goes to zero, see for example ( Conti, 
2003). 
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L(13,2) L(13,2,13) 




Figure 1: Scheme of fields in a structure that realizes the translation bound 



Rank-one path through given field ranges. The described procedure can be iteratively applied to 
find an optimal laminate that links fields a in the given ranges. The process is continued until all materials 
are joined together and the average field of the whole structure is equal to the required value eo- In the 
eigenvalues plane, the process of formation of an optimal structure is represented by a graph that joints 
the permitted values e, by either horizontal or vertical intervals, and the last interval must pass through 
the given point eo. The volume fractions that determine the positions of the fields within intervals are not 
shown in this graph. They must be computed separately, and the constraints on volume fractions determine 
the domain of applicability of the graph. Usually, the graph (and corresponding optimal microstructures) is 
not unique. 

The next two examples demonstrate the earlier obtained optimal three-phase laminates which are parts of 
the whole picture presented here. The structures were originally obtained without the assumption ks = oo, 
but we keep it here. 

Example 2: Optimal structures for the Translation bounds Following (Albin et al. (2007)), 
we construct anisotropic optimal structures that correspond the following sufficient conditions for the fields 
in materials 

ei = [Xai, (1 - A)ai], e 2 = [a 2 ,a 2 ], e 3 = [0, 0] (38) 

where A is an arbitrary real parameter, and oi, a 2 are constants. Field e\ in the first material is un- 
derdetermined, and the field in second material is proportional to a identity matrix. These relations are 
derived from Translation sufficient condition (Albin at al, 2007a). We show here the rank-one path (laminate 
microstructure) that realizes the condition. 

The permitted by translation bounds fields in materials one and three are in rank-one contact if we take 
A = or A = 1 in (|38[) . Joining by laminates certain amounts of materials one (where we take A = 0) and 
three, we obtain laminate (interval Aio) and call it £(13). Fixing the relative volume fractions of materials 
in the laminate, we arrive at the point [0, ei3] = [0,02] that is compatible with material two, field [a 2 ,a 2 ]. 
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Then a second-rank laminate 1/(13, 2) is formed by adding this material in an orthogonal layer, its (average) 
field is represented by a point ei2,3 of the interval [ei3,e2]. Then, we join materials one and three (where 
we take A = 1, see (|38p ) in an orthogonal laminate, choose the relative volume fractions so that the average 
field in the laminate is compatible with the filed ei2,3, and laminate it with the structure 1/(13, 2) obtaining 
structure 1/(13, 2, 13). Assume that all amounts of materials two and three are used to build the structure, 
but there remains an extra amount of material one. To add this material, we choose an appropriate value of 
A in (|38|l . laminating the structure £(13, 2, 13) with the first material and forming structure L(13, 2, 13, 1), 
then laminating it again in an orthogonal direction, obtaining structure 1/(13, 2, 13, 1, 1). The parameters of 
the last two laminations must be chosen to reach a prescribed average field. 

The limits on relative volume fractions or on the attainable by this construction average field can be 
transformed into range of applicability, see (Albin et al. (2007)). It turns out that volume fraction mi must 
be above a threshold g m {r), mi > g m (r) and the rate of anisotropy r must be large enough, r > g r , too. 

Example 3: Isotropic optimal three-material structures Similar technique was used in (Cherkaev, 
2009) to determine optimal structures of isotropic composites (r = 1) beyond the applicability of translation 
bound, mi < (? m (l). The optimal isotropic structures as similar to the structures shown below in Figure^ 
regions A, B, Dl). It is shown that one of three types of laminates are optimal, depending on the range of 
volume fraction mi. There are two critical values < moo < mo < 1 that determine the topology of optimal 
structures: If mi £ [mo, mi], the described above structures 1/(13, 2, 13, 1, 1) are optimal, if mi € [moo, mo], 
the optimal structure degenerate into £(13, 2, 13), and if mi £ [0, moo], the optimal structure becomes 
1/(123,2,123). The regions of applicability of both structures match the regions of applicability of the 
bounds. 

Remark 3.4 Below in Section \6.SX we suggest different than L(123, 2, 123) optimal structures for the case 
mi < moo, that are more convenient for description of anisotropic composites. This variation is possible 
because the rank-one path between fields in optimal ranges is nonunique. 

4 Results 

In this section, we summarize the results. Their derivation is shown in the two following sections. 

Bounds for the energy and minimizers. The analysis in Section [S] demonstrates that the optimal 
bound for energy is a multifaceted surface: it is expressed by different analytic formulas in different domains. 
These domains are conveniently presented in the parameters plane r, mi , Figure [2] where the parameters 
were set as ki = 1, = 3, m,2 = .5 The dependence on m,2 is not shown on the figures. Variation of m,2 leads 
to variation of the shape of the regions above but does not change their topology. It turns out that there 
are five nontrivial cases (Regions A, B, C, D, and E), Figure[2] depending the values of t op t (ki, fe, mi, m,2) 
and r. Region D corresponds to the known translation bound (Milton and Kohn, 1988), optimal structures 
in region Dl have been found in (Albin et al. (2007)), the isotropic structures were found in (Gibiansky 
and Sigmund, 2000). The upper interval r = 1 correspond to isotropy of K*, the isotropic bound have been 
derived in (Nesi 1995), and optimal structures have been found in (Cherkaev, 2009). Other bounds and 
structures are new. 

Boundaries between regions in Figure[2]are calculated, the expressions are shown in Table[T] The division 
between Al and A2 and between Dl and D2 are based on structural attainability, the bounds are given by 
the same expressions. 

Table [2] summarized the expressions for optimal energy bound in different regions, expressions for G- 
closure boundary, and optimal laminates that realize the boundary or best approximate it. Tableland 
Figure show the details of optimal parameters and minimizers in different regions. 
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Figure 2: Regions A-E of multifaceted boundary for the energy in an anisotropic field in the r,m 
plane. 



Tabic 1: Boundaries between regions 



Boundary 
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A2B 


BD 
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BC 


D1D2 


D2E 


Formula 
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(En 
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Table 2: Formulas for energy bounds and G-closure, and type of optimal structure or the best approximate 



Region 


Energy 


G-closure 


Structure 


Exact? 


Al 




(El 


L(13,2,13,2) 


yes 


A2 


fa} 


Cm} 


L(123,2) 


yes 


B 


m 


(52) 


L(13,2,13) 


yes 


C 


m 


(El 


L(13,2) 


yes 


Dl 




(Hi 


L(13,2,13,l,l) 


yes 


D2 


dm) 




L(13,2,l) 


not known 


E 


m 




L(13,2,l) 


no 
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Table 3: Character of minimizcrs 



Region 


"Opt 


d{x) 


in fii 


d(x) in tt 2 


A 


t = k 2 




- d*i2 = Sn 


d 2 i + (^22 ^ ^21 




B 


t € (k u k 2 ) 


d?H 


- d 12 = 5 n 


tfej = 




C 


t 6 (ki,k 2 ) 


dn = 


= Sn, di2 = 


d 2 i = const > 0, c?22 


= 


Dl 


t = ki 




- d 12 < Si! 


d 21 =0 d 22 = 




D2 


t = ki 




- d 12 < Sn 


d 2 i = d 22 = 




E 


te (o,ki) 


dn = 


- const > 0, d\ 2 = 


d 2 \ = const > 0, d 22 


= 




Figure 3: The eigenvalues of the gradient field- minimizers, according to the bounds. The equilib 
rium condition is not assumed. 
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L( 123 ,2) Region A2 L(13,2) Region C 




L(13,2,13,2) Region Al L(13,2,13) Region B L(13,2,13, 1,1) Region Dl 




L(13,2,l) Regions D2, E 



Figure 4: Cartoon of optimal structures in regions A -Dl and the presumed optimal structure ii 
region E. 



Three special points In the map of the regions, notice two points where several region meet. The four 
regions Al, A2, B, and C meet in the point PI: 

fci(l-?n 2 ) , . 

r = m 2 , mi = — (39) 

The four regions B, C, D, and E meet at the point P2 

fei(l-TOa) 

r = m 2 , mi = — ^ : — -. (40) 

fei + k 2 

At the line (PI, P2) the field E 2 is constant and proportional to the identity matrix, and below this line the 
proportionality is lost. 

The three regions A2, C and E meet at the point P3: 

fci(l-m 2 ) , . 

r = 0, m 1 = — K — '-. (41) 

At this point, the L(12,3) structure has the same conductivity in x\ direction, as a simple laminate, see 
(Cherkaev and Gibiansky, 1996). 

G-closure boundaries. The results for the G-closure boundary follow from the energy bounds, the 
expression for G-closure boundary are summarized in Table [21 In FigureJS] the eigenvalues k*i(r), k*%(r) 
are shown that form the G-closure boundary for the values k\ = 1, k 2 = 2, mi = 0.15, m 2 = 0.5. Notice, 
that decreasing values of r correspond to regimes D, B, C, and A (see Figure [2}. The horizontal line in the 
graph represent region C, because both eigenvalues are constant independent of r. 

Figure|B]shows the G-closure boundary, parameter r is excluded. The graphs show the result of comparing 
the bounds obtained here with the translation and harmonic mean bounds for two values of mi . 



5 Derivation of the bounds 

Solving minimization problem (|26p for fixed t, we notice that the algebraic form of Y(r,t) depends on 
sign(fci — t) and sign(fc 2 — 2). There are five cases to analyze: t < k\, t = ki (translation bound), fci < t < 
fc 2 , t = k2, and t > fc 2 . After the solution Y(r,t) of (|26[) is obtained, the bounds are derived by maximizing 
the bound with respect to t, and G-closure bound follows. 
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Figure 5: The eigenvalues of effective tensor A* at the G-closure boundary in dependence of 



5.1 Intermediate value t e (ki,k 2 ) 

Study of minimization problem (I20p If ki < t < fe, V = Vy + V% is computed using (J2S1), ((311)) 

V = 2mifci(Si 2 i + S12) + m 2 (fe + t){S 2 21 + S 2 22 ) + m 2 {k 2 - t)(D 2 21 + D 2 22 ) (42) 
We compute Y(t,r) from the relations (see (|27[) , l|33[l) 



y (t, r) = min V 

Sll, S12, 521, ^22, ^21, -°22 

subject to: 

1 + r 

miSn + m 2 S 2 i — S01 — — m 1 S\ 2 +m 2 S 22 = Q, (43) 
v2 

1 — r 

miDu + miD-ix = D i = — m 1 D 12 + m 2 D 22 = 0, (44) 

V 2 

D'L + D\ 2 < + S 2 21 

(45) 

The problem for Y is a standard finite-dimensional constrained optimization problem for the averages Dij 
and Sij. The analysis of optimality Karush - Kuhn -Tucker conditions leads to the following optimal values: 

fti = S 22 — D 22 = 



Depending on whether or not the last constraint (|33[) is satisfied as an equality or as a strict inequality for 
optimal values of Dij and Sij. We recognize two cases: 
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k1=1,k2=2,m1=0.05,m2=0.5 




23456 4 5 6 7 



Figure 6: The G-closure boundary and its estimates (harmonic and translation bounds), ki = 
l,fc2 = 2, m,2 = 0.5 Left: mi = 0.15 (Regions D, B, C (sharp bend point), A). Right: mi = 0.0. 
(Region A) 
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1. (Case B): If mi Su > -Doi then in the optimal values are: Dn = and D21 = 0. 

mi 

2. (Case C): If < mi-Dn < D01 even when flu reaches its maximum flu = S11, that is equivalent to 

Du = Sn < , then 

r»2 



D21 = > and Dn = Su 

m 2 



(46) 



In this case, the d-minimizer is constant di(x) = -D11 = Su, d^ix) — 0, a.e. x G Qi. 

Now we compute minimum Y(r,t) of V and optimal fields Dij and Sij for these two cases as functions of t 
and given parameters r, mi, and ki. 

Case B If D21 = holds, ([42]) becomes: 

V = 2mifa5n + m 2 (k 2 + t)Slt 

The minimum of V subject to the constraint (|43[) can be found using the standard Lagrange multiplier 
procedure: 

,2 _ (1 + r) 2 



Y B (t,r) = H SZ 



-H 



and the optimal Sn are: 



where 



Si 



l + r 



1 2^2fci 



Ho, S21 



l + r 



V2(k 2 + t) 



H 



(47) 
(48) 



Ho = 



mi m.2 
2ki k 2 + 1 



and point-wise minimizers satisfy 

di — D21 — 0, d2 — D22 — 0, S2 = S22 = 0, si = S21 , a.e. in 0,2 



d\ + d 2 — S11, si = Su, S2 = S12 = 0, a.e. in 51i 



which implies that: 

Tr(e) 

Bound Bb can be found by solving: 
There are two critical values of t: 

icr 1 



cr2 



constant, det(e) = a.e. in Qi, 
e = al, a.e. in 0, 2 



= max (Y B (t,r) - 2tr) 

*e(fci,fc 2 ) 



r mife + 2r m2fci — fci^/rm2(l + r) 
r mi 

r mife + 2r m2fci + fci^/rmjfl + r) 
r mi 



(49) 



(50) 



and it can be checked that maximum is achieved at t crl , therefore t opt = t crl . Substitute it back into Bb 
and we get energy bound: 



fci(l + r-2^ , , 

B B - h r k 2 

Imi 



(51) 
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G-closure boundary After the expression for Bb is obtained, a parametric representation for k*x and 
k,2(r) at the boundary of G-closure are found using (jXTJ) : 

k*i = K(r) k*2 = K[-) (52) 



where 

rst \ fci (1 — Jrm 2 ) (1 + t — 2Jrm 2 ) 
K(r) = — + rk 2 

mi 

Eigenvalues fc*i(r) and fc*2(f) are linked at the boundary of G-closure through parametric equation 

112 

+ ■ 



fe*i — t C ri(r) fc* 2 — t crl (r) Ho — 2t crl {r) 
We show below in Section [6] that bound Bb is achievable by laminate £(13, 2, 13). 

Case C If D21 > 0, V is expressed as: 

V = 2mifciS?i + m 2 (fc 2 + t)Sfi + m 2 (£;2 - t)Dj»i (53) 



substitute D21 given by (|46[) into V and also express S21 in terms of Six using (I43[) we get a quadratic 
function of a single variable Sxx' 

, r „ , c 2 . /, ,x (S01 — mi5n) 2 . . (D01 — miSn) 2 

V = 2mifei5n + (fc 2 + t)- — + (fc 2 — t)- — 

m 2 U12 

Differentiate it with respect to Sxx and we can find the optimal Sxx and Y (t, r) 
„ S x(k 2 + t) D 01 (k2-t) (k 2 + rt) 

oil ~ = 1 = = 1= — (04) 

2k 2k ^/2k 

—mxr 2 t 2 + 2rn 2 (ki — k) rt + m 2 fcifc2 + kk2r 2 

Y c = 2 i 4 

2m 2 fc 

where 

k — mxka + m 2 ki (55) 
Differentiating Yc(r,t) — 2rt with respect of t, we find the optimal value t opt : 



"opt 



(kx - fc) 



nt 2 



rmx 

and the bound 



(56) 



(1 - m 2 ) 2 kx + mim 2 fc2 , k 2 2 

B c = h r (57) 

2mi 2m2 



G-closure boundary According to (JTTJ), we have 



, (l-ra 2 ) fa+mim 2 fe , k 2 . . 

= ' , K*2 = (58) 

mi 7712 

The whole region corresponds to a single point at G-closure boundary. 

Remark 5.1 In optimal two-material mixtures, the optimal two-material structure is a simple laminate 
when r is less than a threshold value. Here, in optimal three-material mixtures, the 1/(13, 2) structure is 
analogous to a laminate, and it also possesses a constant fields in each subdomain Oi . 
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Boundaries of region B and C The applicability of bounds is described through systems of inequal- 
ities of r in terms of conductivities ki and volume fractions rrn. Part of inequalities is obtained by solving 
the constraint fa < t opt < fc 2 . The rest is derived from the requirement 

An „ _ An ^ „ ^ 

< An m case B or > bn m case C. 

mi mi 



The results are presented below. 
Region B 



4>ab < r < ipBD, r > rri2 
fa 

Jk + \Jk 2 — m 2 k 2 
2fa 

+ ^/a 2 - 4m 2 kf 



^ab = m 2 f fcl ] (59) 



/ 2fa 

V'sn = w 2 = 1 , a = k + 2m 2 fa (60) 



Region C 



4>AC < r < IpCD, V < m 2 

m 2 (fa - fc) m 2 (fa - k) 

IpAC = ; , iPce = ; • (61) 

mi K2 mi fei 

The analysis of constraints (|26|l tells that bound Be is valid, if 

£>ii=5u, miDu - Dot < 0. (62) 



otherwise, bound Bb is effective. To work out this inequality, we substitute expression (|54[) for optimal 

1 — r 

-Dii (t, Dqi , Soi) into (1621) . account for relationship _D i = Soi, and obtain condition 

1 + r 



r (mit op t(r) + k^— fam 2 < 



of realizability of region C. Substituting t opt (see (|56p ) into the above inequality and simplifying, we bring 
this condition to the form: r — ?n 2 < in Region C. The boundary between regions B and C corresponds 
to the equality 

r - m 2 = (63) 

5.2 Case t = k 2 . 

Region A: We follow the same procedure as in the previous section. First we find V = Vi + V2 by 
evaluating Vi, V2 in this case (see (l30l) , (p5l) 1: 



V = 2fa ? ni(Si 2 i + S 2 2 ) + 2fami(S 2 2 1 + S 2 2 2 ) 
V is independent of Dij, D 2 j. We minimize ^subject to the only constraints 

miSij + m 2 S2j = Soj, j = 1,2 

The constraints 
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are satisfied by point-wise minimizers dj,j = 1,2 in = 1,2, such that: 

d\ + d 2 < s i + s 2i an d s j = *Saj> a - e m ^2, 

d\ + d 2 = si + si, and Sj = Sij, a.e. in fix, (64) 

These minimizers exist, because Doj < Soj for r G [0, 1]. 

The minimum of V is found using the standard Lagrange multiplier procedure: 

Sij = ^-S j, S 2 j = ^-Soj, J = 1,2 (65) 



2fci J ' * J 2k 



H 2 (l + r) H 2 (l + r) 

Oil — 7= , D21 = 7= , <5l2 = 022 = U 

2^2fci 2y/2k 2 

where 

_ (mi m 2 \ 1 
H2 -\2ki + 2k 2 ) ■ 
We conclude that Su and S21 satisfy the relation 

Sn _ hi_ 

S 2 i ki 

Besides, ()64|) and (|65|l implies that point-wise minimizer in f^2 satisfies: 



(66) 



Tr(e) — constant, a.e. in ft,, i = 1, 2, 

det(e) = 0, a.e. infix, (67) 



1 + r 

Accounting for the value of S01, S01 = — 1=-, S02 = 0; we compute 

v2 

Y A (r) = g a fig 1= H2(1 2 +r)2 (68) 

Bound B \ is: 

B. = ^±^-2^ (69) 



Remark 5.2 VFe notice that when t G (A;2,oo), expression for V is the same as in the case where t — k 2 
(although the pointwise dj, in fi 2 is different, see Section \3.£\ for details). Therefore, expression for Y is the 
same as Ya (see (686 ) and it is independent oft. Maximizing Y — 2tr with respect to t to find the bound, we 
find that t opt = k 2 , because Y A — 2tr linearly decreases with t. In conclusion, if it is assumed that t £ (k 2 , 00), 
the optimal value is t — k 2 and the bound is as in Ba. 

G-closure boundary Find fc*i(r) and k t2 (r) again using (|52|) where 

fc.i = K A (r) = hr + 1)H 2 - rk 2 , fc, 2 = K A (l/r). 

Excluding r from these two equations, we find the equation describing the boundary. One can check that 
fcti(r) and fc* 2 (r) are bounded as 

2 11 

H 2 - 2k 2 ~ k7i -k 2 + k, 2 - k 2 ( ' 

This relation is similar to Translation bound (region D). If the composite is isotopic, the bounds coincide 
with the one found in (Nesi, 1995) 
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Region of applicability The region $a where bound Wa is effective is bounded by inequalities mi > 0, 
< r < 1 and by curves 4>ab and (j>AC- 

r < 4>ab, m 2 < r < 1, 
r < cj>AC < r < m 2 

In coordinates mi,r, the region is shaped as a curvilinear pentagon, It corresponds to small values of mi 
and all range of r and is side neighboring the region B and C. 

5.3 Case t — k\ ( Translation bound) 

Region D This region is analogous to region B. It corresponds to the conditions t opt = fci and moderate 
anisotropy level that corresponds to optimal values D 2 j = 0. The bound corresponds to the classical Hashin- 
Shtrikman (1983) (for r = 1), and translation (Milton and Kohn, 1988), (for r < 1) bounds. The bound 
can be viewed as a special case of Case B, that corresponds to special value of t. The expression for V and 
averages are as in (|48[) where we put t = ki: 

lk\ fci + ki 

D n = £^-<Sn, D 12 =0, D 2j = (71) 
mi 



where 



i mi m 2 



2ki ki + k 2 
We compute Yd = HiS^ and the bound Bd is: 

Cl + r) 2 

go = Hi l 2 ' -2rkx (72) 

Minimizers s(x) are constant in all domains, Sij(x) — Sij. Minimizer d(x) is zero in fl2 and freely varies 
in Q,\ as long as the inequality d 2 (x) < . Sfj holds, keeping its mean value Dij as in (|71[) . 

G-closure Boundary In Region D, the calculations are similar to case A with k 2 replaced by ki and 
H 2 by Hi, see (|70[) . The effective conductivities are given by (|17p . Parameter r can be excluded, and 
the boundary of G-closure corresponding to the Translation bound is (compare bounds for two-component 
mixture Lurie and Cherkaev, 1982, Tartar, 85): 

+ ir^-ir = rs 2 o,. (73) 



fc*i — fci k t2 — ki Hi — 2ki 

5.4 Small values t G (ki,0\ 

Region E corresponds to t opt < ki. In this case, Vi's are convex functionals of d(x) and the minimizers 
d(x) and s(x) are constant 

s(x)j — Sij, d(x)j = in f2j, if < t < fci 
The minimizers are computed by the same procedure as before: 

Sn = °* H+, Dn = 01 H-., Si 2 — Di 2 = 
mi m2 rr _i mi 777,2 



1 = ^ + ^ ) H -i = ™ + 



ki + 1 k 2 + t fci — t k 2 — t 
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We compute bound Be- 

B e = max (Y E - 2tr), (74) 

t6[0,fci) 

Ye = ±[H + (l + r) 2 + H-(l-rf] 

Optimal value t op t of t in region E is a root of the equation d ^ YE ^ 2tT ^ = which is a fourth-order equation 

of t. 4>d 2 e is found by solving - \ t=k 1 = for r. After t op t(r) is founded and substituted into 

1)74( 1 . one can compute the G-closure boundary using ([17)1 . We are not showing this calculation here. 

Summary: Regions D and E The regions are curved tetragons located in the domain of larger mi. 
Region D is adjoined to B along the curve tpBD, 1)60)) and E is adjoined to C along the curve ipcE, 1)61)) . 
These two regions are divided by the curve 

mib — 2mik\k + ay mi (mi — 1) b , . 

<Pdie = = (75) 

a = k + (mi + 1712) ki, 
b = a 2 — (ki + k2) 2 m\ — 4m2fci 2 , 



and are described as 



$d : 4>db < r, 4>die < r < 1, mi < 1 — m 2 , 
>!>£; : max{0, c(>ce} < r < (/>D2E, mi < 1 — m 2 , 



6 Optimal structures 

Here, we describe structures that realize the bounds. All of them are laminates of a rank, obtained by 
sequentially adding to the existing laminate a new one. The normal of laminates always coincide with one of 
mutually orthogonal axis x\ and X2- The fields in the materials is also codirected with these axes everywhere, 
they are denoted as ey = [ct n ,/3 n ] where the first index i shows the location of the layer, second index j 
shows the material, a n and j3 n show the intensities of the field along the x\ and X2 directions in a laminate 
labeled n, respectively. The relations between parameters a n ,/3 n and s n ,d n are given by (|10f) . The plane of 
eigenvalues (a, /3) is the rotated 45° plain of (si,di). It such laminates, dia = ^22 = identically. The field 
in the third material is always zero, which reads — [0, 0]. 

We prove the optimality of the structure by a straight calculation of the fields in materials in an 
optimal structure using the conditions of the rank-one connection - continuity of the tangent component of 
the fields in laminates and the sufficient optimality conditions found in Section [5J. 



6.1 Intermediate Region B 

The multiscale laminates (L(13, 2, 13)-structure) that connect the fields satisfying the sufficient conditions 
shown in Figure [4l region B is attained by the following steps. 

1. L(13) substructure is formed: materials ki and kz is laminated along X2 direction with relative volume 
fraction of material k\ equaling to pii. The fields in the materials are called en and ei3, and ei3 = 0. 
The average field in the substructure is called eio 

2. L(13, 2) substructure is formed: The obtained £(13) structure is laminated with material &2 (field 622) 
along xi direction with relative volume fraction of material equaling to ^2 to form a second rank 
laminate. The average field in the substructure is called e2o- 
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3. L(13, 2, 13) substructure is formed: The obtained £(13, 2) structure is laminated along X2 direction 
with another laminate which is formed of material ki (field e3i) and kz (field e33 = 0) with laminating 
direction parallel to x\ and relative volume fraction of material k\ equaling to [iz\. The relative volume 
fraction of the second rank laminate is [i^. The average field in the substructure is called e3o- 

The average fields in the described substructures are 

eio = en/in e 20 = eio(l — fJ.2) + e 22 /i 2 , e3o = e 3 i/i3i(l — ^t 4 ) + e 20 /i4- (76) 

The laminate's volume fractions satisfy the geometric constraints: 

/Ltn(l — ^2)^4 +/i3i(l — M4) = mi, H4fJ-2 = m 2 (77) 

We compute the fields eij and volume fractions fiij that satisfy both rank-one conditions and sufficient 
optimality conditions (|49p found in Section [5,21 Thus we show the realizability of the bound and find the 
applicability region of the bound. 

Notice that material k\ is always laminated with material k$ with zero field; using the rank-one connection 
- continuity of the tangential component of the fields in laminates, we see that the first entry of en and the 
second entry of e3i are zeros; such fields satisfy sufficient condition (|34|) . The sufficient optimality condition 
(|49|) requires that the the magnitude of e is constant in Hi, therefore 

eu = [0, p], e 3 i = [A 0] 



where is a constant. Optimality condition (|49[) also requires that field in Material fc 2 is proportional to 
identity matrix: 

e 22 = [ai, M ai=/3r (78) 

The remaining group of continuity conditions concerns the average field do in the substructures. Based 
on the rank one connection condition we have the following: 

eio [2] = e 22 [2] -> /3/in = ft (79) 
e 2() [l] = eso[l] -> a?i/i 2 = / 9/i3i (80) 

The average field of the mixture equaling to external field leads us to: e3o = [1, r], hence: 

1 = QiM2(l — Ms) + ^M3lM5 
(^ii(l-/i 2 ) + /3iM 2 )/i4 =r (81) 

Constraints and parameters of an optimal laminate Solving ([77)) . (T78)) - ((8T]l . we obtain: 

(1 + r) — 2Jm 2 r j r 

P = , ai = pi = </ — 

mi y m 2 

rmi 

M11 = -— — . — , = y/rm 2 

(1 + r) x /rm 2 — 2r m 2 

mi / m 2 

M3i = -JT-. — : „ , ! M2 = 1/ 

(1 + r) — 2- v /rm 2 y r 

A straight calculation confirms that the fields coincide with the fields computed for the bounds. 
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Region of applicability Requiring all volume fractions fall into the interval (0,1) and also enforcing 
< r < 1, we have a system of inequalities of r: 

rmi 

< . ■ = — < 1, 

(1 + r)^/rm 2 — 2r m 2 

mi 

< (1 + r) - 2^Fm=; < ' 

< ,/— < 1, < Vrml< 1 
V r 

Solve the above inequalities and we found that: 

{m2 < r < 1 if mi < 2 (s/rn^ — m 2 ) , 

1 - a s - VI - 2a 5 
m2 < r < , otherwise . 
as 

which is consistent with the region of applicability of the bound in region B. 

Region C Note that if r = m%, then /xa = 1, which implies that the inner layer of composite disappears or 
the composite degenerate into T structure - second-rank laminate L(13, 2) Figure 21 region C, that matches 
the bound in region C. This structure plays the same role as laminate in two-phase problem. There, an 
optimal structure degenerates into laminate if r is small enough. 



6.2 Region A 

Region A%: L(13,2,13,2)-structure: Laminates L(13,2,13,2) whose field inside each phase match the 
sufficient conditions shown in Figure|4]region A are obtained by adding a layer of material fc 2 (field e4 2 ) along 
direction x\ to the laminate L(13, 2, 13) described in previous section, with the relative volume fraction of 
1/(13, 2, 13) defined as The average field eo in the structure is again eo = [1, r\. The average fields in the 
substructures are as in (|76[) and additionally we have e4o = e-iofJ-5 + ^42(1 — fj.5). Continuity conditions for 
average field e;o in the substructures are 

eio[2] = e 22 [2], e 20 [l] = e 3 i[L> 3 i, e 30 [2] = e 42 [2] 

Average field e^o — [1, r] leads to: 

1 = a 2 (l-/i S )+/3/i3iM5 (82) 
The volume fractions (relative and absolute) are related as: 

mi — /-tn(l — ^2)a*4/J5 + M3i(l — /J4)a*5, rn-2 — A*2M4A*5 + (1 — Ms) (83) 
The optimality conditions (|66|l . (|67|l yield to relations 

en = [0, P], e 3 i = [/3, 0], e 22 = [«i, Pi], e 42 = [a 2 , r], 
oci + Pi = ct2 + r kiP = fc 2 (ai + Pi) 

where 0:1,02, P, Pi are some constants. The established relations allow for solving for the unknown param- 
eters of the structure - the constants ai,a 2 , P, Pi and 

Remark 6.1 Here, the number of unknowns is bigger than the number of constants. To handle the uncer- 
tainty, we additionally assume that e 22 is proportional to identity matrix, i. e. e 22 [1] = e 22 [2] or 

ai=Pi, (84) 

which significantly simplifies the calculations but might restrict the domain of applicability. 

L(13,2,13,2,2)-structure is optimal in the isotropic case (r = 1) and the field inside material fc 2 in the 
core part is proportional to identity matrix (Cherkaev, 2009). The suggested here L(13,2,13,2)-structure is a 
degeneration of the L(13,2,13,2,2)-structure, so we keep the assumption on the field inside core material fc 2 . 
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Calculation of the constants Solving equations ([82l) - (f84)l . we find the volume fractions and fields 
inside each material: 



where: 



fci a 2 fcf 2a 2 fc 



Mn = 777—, A*31 = — : , = 

2k r(l + r)ai 

= 7—7 — --7T, Ms ~ 



fci(r + l)' ^° (2fc-fci(l + r)) 2 ' 

fe 2 (r + l) rn in fca(r + l) M 

en = — i-= i- [0, 1] , e 3 i = — [1, 0] , 

fe « 

fci(r + l) fci(r+ 1) 

e 2 2 = -~ [1,1], e 4 2 = = 1, 0] + r[-l,l 

2k k 



ai = — 5rm 2 fci 2 — 4r k\m,\ki + Ark 2 + (1 + r — m 2 )fci 2 
a 2 = k\mi(rni — 1) + m\ki(m?. + r) 

Notice that fractions (U 2 and ^131 vanish simultaneously. Such degeneration brings the L(13, 2, 13, 2) structure 
into L(13, 2)-structure 



Region of applicability All the volume fractions have to fall in the the interval (0,1), therefore we 
have the following inequalities: 

g (m 2 fc — fcim 2 + rmifc 2 )fci 2 ^ 
r fc 2 a 2 

2(m 2 fc — fcim 2 + rm 1 k 2 )k 1 ^ 
r a 2 

2k r(l + r)a 2 

0<l— r tt < 1, 0< V ' - < 1 

fcl(r + l) (-rfci - fci +2fc) 2 

The above system of inequalities have solutions: 

*1>A1A2 < r < IpAlB, (85) 

where 

(fci — fc)m 2 m 2 fc 2 

?/>AlA 2 = 7 , tpAlB — — , = 

mi/C2 2fc 2 -m 2 fc 2 + 2fcy / fc 2 -m 2 fc 2 

only if fe, m,i, i = 1, 2 satisfy: 

fci(l - m 2 ) 
mi < 2fc 2 ' 

At the boundary of the applicability domain where r — ipAiB we compute /j, 5 = 1 and conclude that 
the structure degenerates into L(13, 2, 13)-laminate. At the other boundary, when r — ipAiA2 we have 
fj, 2 = 0,^31 = which means that the composite degenerates into £(13, 2)-structure. 

Region A2: L(123,2)-structure In the region A2, the optimal structure is second-rank laminate 
L(123, 2), see figure Region A2. The structure contains the following fields: In the inner layer, ei = {at, 0) 
in fii, ei2 = (c*2,0) in fi 2 i, and es = (0,0) in SI3 , and the in outer layer, e 22 = (022, /S22) in f2 22 . Here Q21 
and Q.22 are the subdivisions of S7 2 , fi 2 = Q. 2 i U $1 22 . 
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The optimality condition requires that a.2 — 0122 + P22 and the compatibility requires that k±ai = k 2 ot2- 
The computation of the constants is similar to the previous case. They are: 

k 2 (l + r) fci(l + r) 

Qi = — a 2 - 



a 22 =r P; 



k k 
fci(l + r) 



22 



k 

The region is bounded by the lines, see Figure [2] 

r = 0. r = 1, mi = 0, t = i\>a\A2, r = ipAC 
When r — > 0, the optimal L(123,2)-structure degenerates into laminate. 



6.3 Region D 



A part Dl of region D is attainable by laminates L(13,2,13,l,l), as it is shown in (Albin et al. (2007)) by 
the method similar to the presented above. The most anisotropic structure of this type is L(13,2,l) with the 
parameters: 



fcl(l — 7712 — rrii 

P= ; '-, en 



mite 



e2 = 



r(ki + A;2)(fci + /c2 — k — miki) 



m2fefa 

r(ki + /c2 — k — m\k\) r(ki + k2 — k — mifci) 



e3i = 



m2te m 2 k 2 

r(l-m 1 )(fci+fc 2 ) 2 r(fci + 2fc 2 ) 



rri2k\k2 



A-2 



and pnii is the volume fraction of material k\ that is rank-one connected with material kz- L(13,2,l) is 
optimal when 



r = 1pDlD2, 
1pDlD2 = 



m,2k 1 k2 



(ki + k2)mi((ki + fe)(l — mi) — 3m2fci) + rri2k\{2k\ + k 2 — 2m,2fci) 
and region Dl is described as 

fci(l - m 2 ) 



mi > 



fcl + k 2 



^>D\D2 < r < 1, r > 'i/'Bo- 



(86) 



(87) 



In Region D2 where r < tpDiD2, we did not find optimal structures. We presume that the translation bound 
is rough in that region. Our guess is to preserve structures L(13,2,l) as the best ones (see discussion in the 
next Section). 



6.4 Attainability of bounds in D and E regions 

The bound (|74[1 in region E is not attainable. Indeed, it suggests that the local fields are constant in each 
material. These constant fields, however, cannot be joined together with a structure. Indeed, an arbitrary 
mixture of fields Ei and E2 cannot be in connection with E3 = 0. Indeed, the determinant of the difference 
of two connected fields must be zero. But one can check by direct calculation that determinant of any convex 
combination of cE\ + (1 — c)Z?2 from (174[) . (|74|l is not zero. The bound cannot be exact: the field in $^3 can 
neighbor neither fields in Q,\ and £^2 nor their mixture. 



2G 



Remark 6.2 An exact bound would require more restrictive constraint than det(e(a;)) > used here. 
These new constraints should depend on volume fractions of materials in the composite or on the mean 
field in it. At present, the needed constraints are not established and the obtained bound is rough (not 
attainable) although very close, see Section\6.4\ 



The part of region D (region D2) close to E is probably not attainable as well. As we described earlier, 
a larger part Dl of region D is attainable by structures are L(13, 2, 13, 1,1), Figure[T] (Albin et al. (2007)). 
In Dl, values of mi and r are larger than the corresponding thresholds. Region Dl adjoins region B along 
curve i/jbd, see figure [21 The most anisotropic structures of this sort are L(13, 2, 1) (Albin et al. (2007)), 
they correspond to minimal values r = 4>d 1 d 2 computed in (Albin et al. (2007)). 



(fei + k-2)m\{{k\ + fe)(l — mi) — 3m2fci) + mik\{2k\ + k 2 — 2ra2fci) 

We have not found optimal structures in the complementary region D2 that is between regions Dl and E 
and conjecture that the bound probably is not exact there. In this region t opt = ki but the bound is not 
attainable. Probably, a better bound is needed for this region and for region E. 



Arguments for presumptive structures We try to guess the best structures in D2- and E-regions. 
The bound is not exact, and we would not attempt to use the method that was exploited for the other 
regions: there is no chance to find the structure that exactly realizes the bound. Instead, we conjecture the 
type of laminate structures based on their asymptotic behavior, adjust an inner structural parameter to the 
external field, and compare the result with the bound, finding the gap between them. 

We notice that at the boundary of neighboring region Dl, optimal structures are of the type 1/(13, 2, 13, 1) 
degenerate on the boundary into the structures 1/(13, 2, 1), Figure[2] In another neighboring region C, optimal 
structures £(13, 2) can be viewed as the result of degeneration of L(13,2, 1), when the volume fraction of 
the exterior layer goes to zero. 

Also, laminates £(12) are optimal two-material (7713 — > 0) structures that correspond to anisotropic 
loading r < rt r < 1, (Lurie and Cherkaev, 1982), these laminates are a degeneration of laminates L(13, 2, 1) 
as mz — > 0. Finally, the optimal structures for the limits when r — > where conjuncted in (Cherkaev 
and Gibiansky, 1988) to be L(13,2, 1), because this structure, being different from simple laminates, has a 
conductivity in x\ direction equal to harmonic mean while the conductivity in orthogonal direction is finite 
(smaller than the arithmetic mean). 

Based on these observation, we presume that structures 1/(13, 2, 1) with the proper distribution of ki 
between the layers, stay optimal in the whole region D2 U E. In those structures, the fields in the second 
and third materials are constant everywhere, as the bound predicts. However, the field in the first material 
takes two different values in different layers; this contradicts the assumption of the bound but makes the 
structure compatible: The field en in the inner layer is rank-one connected with E3. 



Numerical results The numerical experiments are performed to see how well the suggested bounds 
approximate an optimal bound. In all the numerical experiments, conductivities ki of each material is fixed. 
So is the volume fraction rrn of each material. And relative difference between the bound Be and energy 
^1,(13, 2, 1) of L(13,2,l) structures (Figured region E) as follows 



SW re 



min^g^i] W Lil3:2 ,i) (a) - B E 
Be 



are calculated for r £ [0,ro], where ro is the threshold value where the t opt in case E becomes k%. Energy 
Wi(i3,2,i) depends on one parameter a - the relative amount of material one used in the inner layer. We 
choose value a op ± of a to minimize the energy stored in the structure. a opt changes with respect to anisotropy 
level r of external field hence W L ( i:ji2 ,i) ■ The results of one numerical experiment are showed in Figure [7] 
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and the parameters are: mi = 0.2, m,2 = 0.5, m.3 = 0.3, fei = 1, fca = 3. As we can see, the relative differences 
are rather small, are in order of f0 -4 and even in order of 10 -7 as r is very close to value of 0. This is true 
for all fixed mi values which fall inside region E in r — mi plane and 8W re i also changes in the same way 
with respect to r for each fixed mi. 

_4 {Relative difference between W „„ „ and B c ) for t<1 vs r 

x 10 L(u,^,i) t 




ure 7: 6W re i relative gap between the energy of the bounds and guessed structure in region E 
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